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In the variational framework, we study the electronic energy spectrum of massless Dirac fermions 
of graphene subjected to one-dimensional oscillating magnetic and electrostatic fields centered 
around a constant uniform static magnetic field. We analyze the influence of the lateral peri- 
odic modulations in one direction, created by these oscillating electric and magnetic fields, on Dirac 
like Landau levels depending on amplitudes and periods of the field modulations. We compare our 
theoretical results with those found within the framework of non-degenerate perturbation theory. 
We found that the technique presented here yields energies lower than that obtained by the per- 
turbation calculation, and thus gives more stable solutions for the electronic spectrum of massless 
Dirac fermion subjected to a magnetic field perpendicular to graphene layer under the influence of 
additional periodic potentials. 

PACS numbers: 73.61.Wp,73.20.-r,72.80.Rj 

Since the discovery of graphene by Novoselov et. al^, there have appeared many theoretical studies on graphene and 
graphene based nanostructures such as graphene quantum dots and graphene nanoribbons. There are also few reviews 
which have been devoted to these structures^. In most of these works, theoretical studies on the electronic energy 
spectrum of these structures have focused, due to their linear dispersion relation near K points in the Brillouin zone^, 
on the continuum version of tight binding Hamiltonian, i.e., 2 + 1 Dirac- Weyl equation^, which has been confirmed 
by the experiments^ 

The theoretical considerations of Weiss oscillations in an electrically modulated graphene were discussed by Peeters 
and Matulis 7 . They showed that, within the framework of non-degenerate first order perturbation theory , these 
non-relativistic oscillations are more pronounced in graphene as compared with those found for its 2D non-relativistic 
counterpart^. Very recently, studies on these oscillations in a magnetically modulated graphene have also been 
reported^—. 

On the one hand, from the non-relativistic point of view, after Weiss et al*^, theoretical studies on magnetoresistance 
oscillations in two-dimensional electron gas structures sub jected to a periodic potential have been intensively studied 
for a long time within the frame work of perturbation theory^. It is shown that these oscillations have as a common 
origin the oscillating bandwidth of the modulation-broadened Landau bandsi£. On the other hand, role of one- 
dimensional periodic potentials is crucial itself in graphene. Very recently, it was shown that, by a Muffin- Tin type^ 7 - 
one-dimensional periodic potential, velocity of charge carriers can be controlled in graphene. These type potentials 
allow one to create various types of graphene superlattices and to control the flow of charge carriers in an non-invasive 
way. 

In this paper, we present a variational analysis of the combined effects of electric and magnetic potentials on 
the electronic energy spectrum of a Dirac- Weyl like electron in graphene. Our method is based upon the use of 
solutions of massless Dirac fermions in an uniform magnetic field as trial wave functions in the presence of external 
potentials. Analytical results obtained in this paper indicate that the variational method is more efficient than the 
other previously used methods, i.e., non-degenerate first-order perturbation theory. 



The time independent effective massless Dirac- Weyl Hamiltonian we consider is H = Hq + U (rj_), where 
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is the Hamiltonian for an electron which is minimally coupled to magnetic field B =(0, 0, B + Bq cosKx) with K = 
2ir/ao, through the vector potential A—(0,Bx + (B o /K)sinKx,0). Here, U (r±) is the one dimensional periodic 
electrostatic potential, and is given by U (r±) = U cos (K'x) with K' = 27r/c'.Here, ao and d are the periods of 
magnetic and electrostatic modulations, respectively. In Eq.Q, we have used the Dirac-Pauli representation of Dirac 
matrix a which is written of two by two-block form in terms of Pauli spin matrices <x. Using two component spinor 
(pseudospin) as ^ = (4>* x*) we see that each component of eigenvalue equation HtfS! = Ety satisfies the following 
coupled first order differential equations: 
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where = E/vp. Decoupling them in the absence of external potentials, Eq. ([2]) yields the second-order differential 
equation 
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where <x 3 is the third component of the Pauli matrices. Since H ,p y — 0, we set 
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exp (ik y y) <j)(x). 
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Therefore, Eq. (j3|) reduces to the solution of second order equation for two-component wave function 4> (x) 

2 



eB 

ky + —X 



eB 
lie 



o- 3 



(j)(x) = 0. 



(5) 



It is easy to show that, in pseudospin basis, the Hamiltonian H in two-dimensions has solutions which can easily 
be expressed in terms of the Hermite polynomials by just setting with A = {he/ eB) 
Therefore, the total solution of Eq. @ can be written as 
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and with 



for the eigenvalues 
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which, in fact, corresponds to the eigenenergy E^' ' — ylnjl for the it* band of graphene, in which the quantum 
number n takes on the values n = v + (s + 1) /2. One can follow the same procedure to obtain the associated 
two-component solutions of Eq. (fTJ) for the hole part of spectrum, i.e., the 7r band of graphene corresponding to the 

eigenenergy E^ n ^ — —\/2n/£. In Eqs. ([S][7]), we have defined x — chk y /eB, and rescaled the energy E by dividing 

it by h to obtain E = E jh = E/Hvf, and defined j 2 — eB/hc, £ — TicjeB. Then, for the value s = 1 we obtain 
v = n — 1, and for s = — 1 we obtain v — n. It should be noted that, s = ±1 does not represent "spin up" and 
"spin down", but they describe states on the A(B) sublattice of graphene. In the presence of external potentials, to 
calculate their combined effects onto eigenvalues given by Eq. (J8]), we choose Eq. ((6]) and Eq. (J7J as basis functions 
for our variational procedure. Therefore, in this sense one has to minimize the energy, 
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where U (rj_) = U (r±) /hvp. Analogously, one can follow the same treatment for the negative energy pseudospinors 
to obtain E (7). Finally, performing the related integrals in Eq. ([9]), we find that 
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where the energy E n j; (7) is again rewritten, for convenience, in the units of Hvf, and all the lengths are in the units of 
ooj 7 — 7 a 0: £ — £/(Xq, £q — Iq/clq, k y = k y ao. Here, we have also defined p = clq/c'. Therefore, our rescaled variables 



in Ecj. (|10| measure lengths in the units of a>o and energies in the units of Hvf- Eq. (|10)) , when minimized with respect 
to 7, gives the effect of external one dimensional electrostatic potential together with its magnetic analogue onto the 
Dirac-Weyl like Landau levels given by Eq. ©. To show this, we start from the well-known case, absence of external 

fields, i.e., Bq, U. In this case, minimization of Eq. ([1 
result back into Eq. (TTU|) yields 



with respect to 7 yields 7 = eB/hc — l/l . Replacing this 
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by using first-order perturbation correction. 
= -\-\f2nfl in the absence of external fields. 



It is easy to show that Eq. (fTTj) reduces exactly to those found in Refs. 
Of course, it covers inherently well-known Dirac like Landau levels, E~* 

In FIG. HJa), we plot dimensionless magnetic confinement length £ variation of the dimensionless half bandwidth 
E n o — y/2n/£ for two different values of U = 0.5 and 1 resulting from the numerical minimization of Eq. (|10[) (solid 
lines), together with the results of non-degenerate first-order perturbation calculation (dashed lines) as well as their 
difference AE (lower panel). In order to observe the effect of one dimensional magnetic analogue in FIG. [lib) we 
also plot the same in FIG. (Ha) but for \ — 1 (x = B /U) with p = 0.5, 1, and 2, respectively. From both figures, as 
we decrease £ ,i.e, increase B = B/B (B = hc/eafy Landau levels broaden into mini subbands whose boundaries are 
determined by two asymptotic values of k y , i.e, by k y = and k y = l/2n£. Of these only with n = 1 and n = 2 are 
depicted in FIG. Ufa) for k y = 0. 

It is clear from FIG. [lja) and (b) that, as £ decreases, dimensionless half bandwidth of the first two Landau 
levels begins to highly oscillate, due to the oscillatory nature of Laguerre polynomials, passing through degeneracy 
restoring points, i.e., wherein the flatband condition is fulfilled, and thus Weiss oscillations are suppressed. It should 
also be noted that, from the lower panels, a rapid increase of AE with increase in U is evident, yielding significant 
discrepancies above O.IHvf- Moreover, note that AE diminishes at some values of £ where approximately flatband 
condition is fulfilled. In other words, since diminishing AE indicates that the variational and the perturbational results 
are almost same, we can say that, around these points, i.e., where the flatband condition is fulfilled, perturbational 
results can be safely used, but otherwise they cannot. 

We further plot in FIG the variations of dimensionless variational parameter 7 as a function of dimensionless 
magnetic confinement length £ for the first two Landau levels with different p values. Also, for comparison, the curve 
for the unperturbed case, i.e., 7 = 1/1, is plotted (dashed bold line). We note two features of these curves. First, 
there exist discrepancies between the unperturbed and perturbed curves in the high magnetic field regime, £ < 1 
where the variational picture becomes more appropriate. Second, again in this region the perturbed curves intersect 
the unperturbed one wherein the flatband condition is fulfilled, as is indicated above. These justify the validity of our 
variational procedure. 

In FIG. [31 E n Q as a function of vB are plotted for the first four Landau levels. With increasing B strong deviations 
from the unperturbed Dirac-Weyl like Landau levels (dashed ones) appear near B ~ 1. To see these effects more 
clearly, variations of half-bandwidths of the corresponding levels, i.e., AE n0 = E n0 — V 2nB, are given as a function 
of y/~B in the lower panel of the same figure. 
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Finally, in FIG. SJa) and (b) we have plotted E n0 as a function of for the four p values, and for the four \ 
values, respectively. It should be noted that the variable p = ao/c' measure whether oscillating fields arc in phase 
or out of phase, and the variable x measures the ratio of amplitudes of fields. One observes very clearly that, by 
changing their phases and/or their ratios of amplitudes, dramatic changes in the spectrum occur when B decreases. 

In conclusion, since our variational method yields energies lower than that obtained by the first-order perturbation 
calculation, it provides an efficient procedure for finding out the effects of external potentials onto the electronic 
spectrum of a graphene electron in an uniform magnetic field. Moreover, we see that, after performing various 
analytical checks on the obtained energy spectra, our results cover the well-known results found in the literature^, 
and in the absence of external fields they reduce to the 2D-Landau levels of massless Dirac fermions in some certain 
values of variational parameter. We use the solutions of massless Dirac fermions in an uniform magnetic field as trial 
wave functions to obtain the low-lying energy spectrum of these particles in the presence of external potentials. For 
Landau bands with high indices in the case of high magnetic fields, adjacent Landau bands may overlap. In this case, 
the variational function should be a linear combination of these sates rather than one of them. In this regime, as 
a consequence of overlapping, even-odd transitions in the Shubnikov-de Haas oscillations^ may be expected in the 
graphene. 
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FIG. 1: (a) E no — y/2n/£ as a function of dimensionless magnetic confinement length 1 — £/ao for the first two Landau levels, 
i.e. for n = 1, 2, in the absence of magnetic modulation for U — 1 and U = 0.5 with k y — 0. While the solid curves represent 
the present theoretical calculations, the dashed ones corresponds calculations using first order non-degenerate perturbation 
theory presented in Ref . 0. In upper panel, while upper curves represent the U = 1, lower curves corresponds to the U = 0.5. 
(b) The same as in (a), but for U — 1 and x = 1- in both figures, the two theoretical curves are compared with each other by 

— — V — PT — — V 

plotting their difference A_B n o = E r j: y — E n ^ y as a function of dimensionless magnetic confinement length £ where E n ~j: y and 

p'Y 

E n k y denote the results of present work and non-degenerate perturbation theory, respectively (lower panels). 
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FIG. 2: The dimensionless variational parameter 7 as a function of dimensionless magnetic confinement length I. While the 
solid lines correspond to energy level with n = 1, the dotted lines correspond to energy level with n = 2, for various sets of p. 
The inset shows the Bq = case for [7=1. 
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FIG. 3: (a) The first four low-lying Landau levels as a function of Vi? with U = 1 for Bo = 0, and for k y = 0. While the dashed 
curves correspond to unperturbed Landau levels, i.e.,V 2nB, the solid ones denotes the effect of electrostatic potential, (upper 
panel). To demonstrate the effect of external modulated electrostatic potential E n o — V2nB is also plotted as a function of 
^fB, in lower Panel, (b) The same as in (a), but for (7 = and x — 1- 




FIG. 4: The same as in FIG. [2] but for various sets of p and \i an d only f° r the first three energy levels. The bold dashed 
straight lines refer to the relevant unperturbed Landau levels 



